THOMAS-FERMI THEORY - Sometimes 
called the 'statistical theory', it was invented by 
L. H. ThomasfTHl and E. FermipF), shortly after 



Schrodinger invented his quantum-mechanical wave 
equation, in order to approximately describe the elec- 
tron density, p(x), x £ R 3 , and the ground state 
energy, E(N) for a large atom or molecule with a 
large number, N, of electrons. Schrodinger's equa- 
tion, which would give the exact density and energy, 
cannot be easily handled when N is large. 

A starting point for the theory is the TF energy 
functional. For a molecule with K nuclei of charges 
Zi > and locations Ri € R 3 (i = 1, K), it is 
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U — Z,Z, It, It; 1 , 

l<i<j<K 

and 7 = (37r 2 ) 2 / 3 . The constraint on p is > 
and J R3 p = N . The functional p — > £(p) is convex. 

The justification for this functional is this: 

• The first term is roughly the minimum quantum- 
mechanical kinetic energy of N electrons needed to 
produce an electron density p. 

• The second term is the attractive interaction of 
the N electrons with the K nuclei, via the Coulomb 
potential V . 

• The third is approximately the electron-electron 
repulsive energy. 

• U is the nuclear-nuclear repulsion and is an 
important constant. 

The TF energy is defined to be 

E TF (N) = m{{£(p) : p e L 5 '\ Jp = N,p>0}, 



i.e., the TF energy and density is obtained by mini- 
mizing £(p) with p G L 5 / 3 (R 3 ) and f p = N. The 
Euler- Lagrange equation, called the Thomas-Fermi 
equation, is 
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where [a]+ = max{0, a), p is some constant {Lagrange 
multiplier) and $ is the TF potential: 



= V(x) 
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\x - y\ 1 p{y) dy. (3) 



The following essential mathematical facts about the 
TF equation were established by E.H. Lieb and B. 



Simon [LS (c/. the review article [EL]) 



1. There is a density pj^ that minimizes £{p) if and 
only if N < Z := Ylj=i This pj^ is unique 
and it satisfies the TF equation (2) for some p > 
0. Every positive solution, p, of (2) is a minimizcr 
of (1) for N = J p. If N > Z then E TF (N) = 
E TF (Z) and any minimizing sequence converges 
weakly in L 5 / 3 (R 3 ) to pY ■ 

2. $(x) > for all x. (This need not be so for the 
real Schrodinger p.) 

3. p = p(N) is a strictly monotonically decreasing 
function of N and p(Z) = (the neutral case), 
p is the chemical potential, namely 
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E (N) is a strictly convex, decreasing function 
of N for N < Z and E TF {N) = E TF (Z) for 
N > Z. If N < Z, pjf has compact support. 

When N = Z, (2) becomes 7p 2 / 3 = <£>. By applying 
the Laplacian A to both sides we obtain 
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which is the form in which the TF equation is usually 
stated (but it is valid only for N = Z). 



An important property of the solution is Teller's 
theorem 



(proved rigorously in [LSQ which im- 
plies that the TF molecule is always unstable, i.e., 
for each N < Z there are K numbers Nj £ (0, Zj) 
with V. Nj = N such that 



E (N) > 
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where E^ om (Nj,Zj) is the TF energy with K 
1 . '/ = Zj and N = Nj. The presence of U in (1) 
is crucial for this result. The inequality is strict. Not 
only does E TF decrease when the nuclei are pulled 
infinitely far apart (which is what (4) says) but any 
dilation of the nuclear coordinates (Rj — > lRj,l > 1) 
will decrease E TF in the neutral case (positivity of 
the pressure) |EL[ , [BL]. This theorem plays an im- 
portant role in the stability of matter. 

An important question concerns the connection be- 
tween E TF (N) and E Cl (N), the ground state energy 
(= infimum of the spectrum) of the Schrodinger op- 
erator, H , it was meant to approximate. 
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which acts on the antisymmetric functions 
A Ar L 2 (R 3 ;C 2 ) (i.e., functions of space and spin). It 
used to be believed that E TF is asymptotically exact 
as TV — > oo but this is not quite right; Z — » oo is also 
needed. Lieb and Simon [LSJ proved that if we fix 
K and Zj/Z and we set Rj 
R° € R 3 , and set N 



lim E 



Z-V 3 R°j, with fixed 
XZ, with < A < 1 then 

TF {XZ)/E Q {XZ) = 1 . (5) 



In particular, a simple change of variables shows that 
El F oln (X,Z) = Z 7 /3£TF m(Ajl) and hence the true 

energy of a large atom is asymptotically proportional 
to Z 7 / 3 . Likewise, there is a well-defined sense in 
which the quantum mechanical density converges to 

TF 
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~ TF (of pi). 



The TF density for an atom located at R = 0, 
which is spherically symmetric, scales as 



p!u>m( x : N 
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Thus, a large atom (i.e., large Z) is smaller than a 
Z = 1 atom by a factor Z^ 1 / 3 in radius. Despite 
this seeming paradox, TF theory gives the correct 
electron density in a real atom — so far as the bulk 
of the electrons is concerned — as Z — > oo 

Another important fact is the large |x| asymptotics 
of pJtom for a neutral atom. As |x| — » oo, 



pItom( x : N 
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independent of Z. Again, this behavior agrees with 
quantum mechanics — on a length scale Z _1//3 , which 
is where the bulk of the electrons are to be found. 

In light of the limit theorem (5), Teller's theorem 
can be understood as saying that as Z — > oo the 
quantum mechanical binding energy of a molecule is 
of lower order in Z than the total ground state energy. 
Thus, Teller's theorem is not a defect of TF theory 
(although it is sometimes interpreted that way) but 
an important statement about the true quantum me- 
chanical situation. 

For finite Z one can show, using the Lieb-Thirring 
inequality 



and the Lieb-Oxford inequality [LO|, 
that J B TF (iV)7with a modified 7, gives a lower bound 
to E®(N). 

Several 'improvements' to Thomas-Fermi theory 
have been proposed, but none have a fundamen- 
tal significance in the sense of being 'exact' in the 
Z — > oo limit. The von Weizsacker correction con- 
sists in adding a term 



(const.) / Vv/p(x)| 2 dx 

to £ (p). This preserves the convexity of £{p) and adds 
(const. )Z 2 to E TF (N) when Z is large. It also has 
the effect that the range of N for which there is a 
minimizing p is extend from [0,Z] to [0,Z + (const.) 
K]. 

Another correction, the Dirac exchange energy, is 
to add 

-(const.) / p(x) i/3 dx 
Jr. 3 

to £{p). This spoils the convexity but not the range 
[0,Z] for which a minimizing p exists cf. [LE] for both 
of these corrections. 
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When a uniform external magnetic field B is 
present, the operator — A in H is replaced by 



|tV + A(x)\' 



B(x) 



with curl A = B and a denoting the Pauli spin ma- 
trices. This leads to a modified TF theory that is 
asymptotically exact as Z — > oo, but the theory de- 
pends on the manner in which B varies with Z . There 
are five distinct regimes and theories: B <C Z 4 / 3 , B ~ 
Z 4/3 jZ 4/3 «B«Z3 ) 5^23 1 »23, These theo- 



ries [LSY1|, [LSY2] are relevant for neutron stars. 
Another class of TF theories with magnetic fields is 
relevant for electrons confined to two-dimensional ge- 
ometries (quantum dots) ]LSY3| . In this case there 
are three regimes. A convenient review is [LSY4]. 

Still another modification of TF theory is its ex- 
tension from a theory of the ground states of atoms 
and molecules (which corresponds to zero tempera- 
ture) to a theory of positive temperature states of 
large systems such as stars (cf. [JM], [ WT| ). 
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